The Kondo effect in bosonic spin liquids 
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In a metal, a magnetic impurity is fully screened by the conduction electrons at low tempera- 
ture. In contrast, impurity moments coupled to spin-1 bulk bosons, such as triplet excitations in 
paramagnets, are only partially screened, even at the bulk quantum critical point. We argue that 
this difference is not due to the quantum statistics of the host particles but instead related to the 
structure of the impurity-host coupling, by demonstrating that frustrated magnets with bosonic 
spinon excitations can display a bosonic version of the Kondo effect. However, the Bose statistics 
of the bulk implies distinct behavior, such as a weak-coupling impurity quantum phase transition, 
and perfect screening for a range of impurity spin values. We discuss implications of our results for 
the compound CS2C11CI4, as well as possible extensions to multicomponent bosonic gases. 



Magnetic impurities can serve as powerful probes of 
the elementary excitations of the host system in which 
they are implanted. The well-studied situation of a spin- 
1/2 impurity antiferromagnetically coupled to electronic 
quasiparticles leads to the Kondo effect, a screening of the 
magnetic moment at low temperature [1]. Subsequently, 
embedding of impurities in strongly interacting quantum 
systems, like high-temperature superconductors, quan- 
tum antiferromagnets, or pseudogap metals has attracted 
a lot of experimental and theoretical work [2-8] . 

Models of impurity moments interacting with collec- 
tive modes of quantum magnets turn out to be par- 
ticularly interesting. Consider a spin system, such as 
the Heisenberg bilayer [9], that can be tuned through a 
continuous quantum phase transition (QPT) between a 
Neel-ordered and a paramagnetic ground state. Close to 
the transition, critical fluctuations can be modelled by a 
standard 0(3) <f> theory. The quantum-disordered phase 
has gapped spin-1 bosonic excitations, a^ a (a = x,y, z), 
which represent the dynamics of the bulk order param- 
eter, <j) a oc (a a + ajj. The magnetic impurity S is cou- 
pled to the order parameter at the impurity site through 
Wimp = lS ■ 0(0). As the fluctuations are gapped, the 
low-energy properties of the impurity are only weakly 
influenced by the bulk, and the impurity susceptibility 
displays Curie behavior, Txi mp (T) = S(S + 1). Surpris- 
ingly, such behavior, albeit with a renormalized Curie 
constant, survives at the host quantum critical point 
(QCP) where the host gap vanishes: an impurity cou- 
pled to a bath of gapless (critical) spin-1 bosons is not 
completely screened [6]. This is in stark contrast to the 
fcrmionic Kondo effect, described by the Hamiltonian 
Wimp = JS ■ 4(0)f (7(7 'C (T ,(0)/2, where f are the 

Pauli matrices and c£(0) creates a conduction electron 
with spin er =|, J, at the impurity site. Physically, the two 
situations are different because the bosonic field 0(0) acts 
as a (classical) fluctuating magnetic field which is unable 
to fully screen the moment, whereas the fermions can 
strongly bind to the impurity to form a quantum singlet. 

The purpose of this paper is to show that this funda- 
mental difference in the impurity behavior is not due to 



the quantum statistics of the bulk particles, but instead 
related to the structure of the impurity-host vertex. We 
shall demonstrate that bosonic host systems, such as frus- 
trated critical magnets with a spin-liquid ground state 
and elementary spinon excitations, can quench an ex- 
tra moment at low temperature, showing a true Kondo 
effect [17]. We introduce a new class of quantum impu- 
rity models, with superficial similarity to the fermionic 
Kondo model, but with bosonic spin- 1/2 bulk particles. 
Remarkably, the bosonic nature of the host implies now 
markedly different properties at both weak and strong 
coupling, as we prove using scaling analysis, solvable toy 
models and a full large- N solution. An experimental mo- 
tivation lays also in the insulating spin system CS2CUCI4, 
in which broad neutron scattering spectra have been at- 
tributed to the presence of spin-1/2 excitations [10, 11]. 

Frustrated magnets. For certain classes of frus- 
trated magnets with half-integer spin per unit cell, a 
4 theory fails to describe the critical properties near 
the order-disorder transition. In the following, we will 
focus on a special kind of spin liquid, in which low- 
energy excitations are deconfined spin-1/2 bosons, which 
were introduced in relation with two classes of models: 
anisotropic triangular magnets with incommensurate cor- 
relations [11-13] (which model the compound CS2CUCI4), 
and frustrated magnets with spin-1/2 per unit cell [14]. 
In both cases, spinons are liberated at the QCP sep- 
arating the paramagnetic from the Neel-ordered phase 
(although the case considered in [12] shows deconfined 
spinons in the whole paramagnetic phase). The effective 
low-energy theory for the bulk QPT takes the form 
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where z a are spinon (or CP 1 ) fields with linear dispersion 
and constrained by |z CT | 2 = 1. Here, gauge fluctua- 
tions are assumed to be gapped and are neglected in (1), 
which is only vindicated [12] for the QPT between in- 
commensurate order and spin liquid in the anisotropic 
triangular lattice (gapless non-compact gauge fluctua- 
tions are important in the critical behavior proposed 
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in [14]). The above constraint implies nevertheless strong 
interactions among the bulk bosons. Importantly, the 
z a are fractional objects in terms of the physical or- 
der parameter: <f> — ^2 acr , z\T aa i z a , /2. This implies 
a large bulk anomalous dimension r/ for the order pa- 
rameter, i.e., the spin susceptibility at criticality follows 
X _1 (w, k) oc (uj 2 — k 2 ) 1 " 1 */ 2 , with t] = d — 1 at Gaussian 
level (note that the standard 0(3) <p A theory has a tiny 
anomalous dimension in d — 2). With the bulk theory 
(1) the simplest SU(2)-invariant impurity action reads: 

Sj = [dr JS -J24(0)^f^(0) + S Bcil . y [S}. (2) 

Note that certain strong-coupling effects, i.e. high-energy 
bulk excitations created locally by the impurity (like non- 
trivial vortex fluctuations), are neglected. 

Weak- coupling analysis. Bare perturbation theory in 
J, neglecting bulk interactions, yields logarithmic diver- 
gencies in d = 2. This suggests a weak-coupling renor- 
malization group (RG) treatment. Importantly and in 
contrast to the fermionic Kondo problem, a potential 
scattering term, S v = J dr{V/A) £ ff 4(0)^(0), is gen- 
erated in the process of renormalization. At two-loop or- 
der, we find the following field-theoretic beta functions: 
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ej - vj ~ \f + 0(f) , (3a) 
ev-\v 2 -\f + 0{f) 7 (3b) 



where J = (2/g)fi~ e j, V = (2/g)[i~ e v, fi is the renor- 
malization scale, and e = d — 2 (for non-interacting bulk 
bosons). Notice that (3(j) does not contain even pow- 
ers of j [18], as the action is invariant under J — ► — J 



and z„ 



furthermore the j term is relevant 
for bosons. The solution of (3) is depicted in Fig. 1: For 
e > it shows a QCP between a local-moment (LM) 
phase, J = v = 0, and a presumed strong-coupling phase 
(to be analyzed below); the RG flow is somewhat simi- 
lar to the fermionic pseudogap Kondo model [7, 8] with 
density of states p{u>) oc \cj\ e . For e < a transition is ob- 
served only in the presence of non-zero bare v (which sta- 
bilizes a bound state, denoted BS), otherwise run-away 
flow to strong coupling obtains. Finally for e = 0, the 
QCP disappears, leaving only flow towards strong cou- 
pling, not unlike the fermionic Kondo problem in metals. 
These results are strikingly different from the behavior 
of impurities coupled to gapless spin-1 bosons [6], where 
the coupling 7 flows to a stable intermediate-coupling 
fixed point, rendering a phase transition and true Kondo 
screening absent. 

Incorporating bulk interactions modifies the scaling di- 
mensions of j and v, as these are determined by the lo- 
cal bulk spin correlations, x( T ) °c l/T d ~ 1+? '; this gives 
e = (d- 3 + r/)/2 in Eqs. (3). Notice that the QCP in the 





FIG. 1: (color online) Impurity RG flow of the bosonic Kondo 
model in dimension d — 2 + e (for free bulk bosons). For e > 0, 
the critical line crosses the v = axis, and a quantum phase 
transition is possible. In the presence of bulk interactions, the 
bulk anomalous dimension dictates the weak-coupling behav- 
ior, so that in d — 2 figure a) [b)] applies to n > 1 [n < 1] . 

incommensurate magnet has 7/ > 1 in d — 2 [11, 12], i.e. 
a boundary quantum phase transition will occur, Fig. la. 
Unfortunately, a consistent RG including the bulk con- 



straint ~^2 a 



= 1 cannot be performed (in contrast to 



the theory of Ref . 6) . Interestingly, two transformations 
allow to map our original model onto a theory where 
both bulk and boundary interactions can be perturba- 
tively controlled: (i) relax the hard constraint on the z a 
in favor of a repulsive u\z a | 4 term for the spinons; (ii) per- 
form a reversed Schricffcr- Wolff transformation as done 
in Ref. 8 for the pseudogap Kondo model. This leads to 
an effective bosonic Anderson model, where both u and 
the hybridization term are marginal in d = 3 which allows 
for an epsilon expansion; details will appear elsewhere. 

Strong- coupling analysis. Let us turn to the strong- 
coupling limit, and ask whether screening can really oc- 
cur at low temperature. Unfortunately, even when reduc- 
ing the problem to a single site, our bosonic Kondo model 
remains complicated. However, we can introduce solvable 
toy models that illustrate the general strong-coupling as- 
pects. In the first example, we replace the z G spinons 
by canonical bosons a CT , described by the Hamiltonian 

Hi = c/(E CT 40 2 + Js ■ Eaa' 4vv/2. This is 

trivially diagonalized, and shows two crucial new aspects 
of the bosonic case which lack in the fermionic Kondo 
model: (i) without bulk interaction (U — 0), the model 
is ill-defined in the sense that no Pauli principle prevents 
the "collapse" of an infinite number of bosons onto the 
impurity; (ii) when U > J/2, screening is possible for 
a range of impurity spin values S < S c = J/(2U — J), 
and underscreening obtains otherwise. This is due to the 
fact that spin- 1/2 bosons can recombine into larger spin 
objects, whose size is bounded due to the local U term. 

Next we verify that these results also apply to the ini- 
tial model based on the z a spinons, as they are not canon- 
ical bosons [cf. the dynamic term in Eq. (1)]. We use a 
second toy model, now based on the z a bosons, where the 
initial 0(4) symmetry is replaced by 0(2) x 0(2) using a 
simplified constraint \z a \ 2 = 1. This is resolved by setting 



3 



z a = exp(i(f> a ), and the one-site bosonic Kondo model 
thus becomes H 2 = -g(d/dO) 2 + (J/2) exp(i0)S _ + h.c, 
where 8 = <f>^ — and S ± = S x ± iS y . H 2 is diag- 
onalized by a unitary transformation H 2 — XJH^U^ 1 , 
with U = exp(i9S z ). For an impurity spin £ = 1/2, 
the ground-state wave function is found to be unique 
and reads = exp(-i0/2) |T) -exp(i0/2) demon- 

strating that the z a spinons can succeed in forming a 
zero-entropy state with a quantum spin. 

Large-N approach. In order to put the previous ar- 
guments in a global perspective, we now develop a non- 
perturbative large- TV approach. The present problem can 
be generalized to a larger symmetry, by choosing S to be 
in the conjugate representation of SU(A) (in order to 
allow for SU(A r ) singlet formation). Eq. (2) now reads: 

Sj= [dT^bUdr-»)b„-^bib a ,4(o)zAo) (4) 

J ° a <j<j> 

where Schwinger bosons b a (a = 1 . . . N) , satisfying the 
constraint b\b a = SN, have been introduced (fi is 
the related Lagrange multiplier). Alternatively, we could 
have taken a full Schwinger boson solution of the bulk 
plus impurity problem, which - after integrating out 
high-energy bulk modes - yields a similar effective theory 
Now, the interaction term in Eq. (4) can be decoupled. 
At N = oo, the problem is controlled by a static sad- 
dle point, with three mean-field parameters: the chem- 
ical potential /i, a bond variable Q = (zj(0)6j.) which 
characterizes magnetic correlations between the bulk and 
the impurity, and a bulk mass term A(x) conjugate to 
zj.(x)z cr (x) which enforces the constraint on the z a . No- 
tice that a space dependence of A(x) is induced by J 
around the impurity. From (4) we obtain the effective 
impurity Green's function of the b a bosons: 



a) 



G 



h (iv n ) =iu n + fi + Q [G z (iv n ) - G z (iO)] 



(5) 



(v n = 2im/P) where the renormalizcd chemical potential 
ft < and the bond parameter Q are determined by: 
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QG z (iv n ) 



G b {iv n ) 



{ -Q/J } 



(6) 



We have also defined the local bulk propagator G z {iv n ) = 
G z (iv n ,0,0), where G- 1 (ii/ n ,x,x / ) = g~\v 2 n - + 
A(x)]<5(x— x'). The spatial dependence, induced by A(x), 
has to be determined self-consistently. Finally, the local 
susceptibility is given by Xioc^) = TY, n G h (iv n ) 2 . 

The form of the propagator (5) allows to predict three 
distinct zero-temperature phases: (i) a local-moment 
phase (LM) with Q = jl = 0; (ii) an underscreened phase 
(US) with Q 7^ 0, jl = 0; (hi) an exactly screened phase 
(ES) with Q 7^ 0, fl 7^ (this nomenclature will become 
clear shortly). The onset of a non-zero Q upon lowering 
T defines the Kondo temperature Tk- To determine the 
actual phase diagram as a function of the parameters J, 
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FIG. 2: Schematic phase diagrams at N = oo for the Kondo 
problem with bosonic spinons in dimension d = 2 + e. 

S, g, we need to solve the full set of saddle-point equa- 
tions. As a first step, the spatial dependence of A(x) will 
be ignored (A(x) — > A), which is justified in the regime 
where Q is small (as discussed below). The local bulk 
propagator is then given by: 



G z (iv n ) 



I 



d d k 



G z (i0) 



X + k 2 
AW 2 +\\ d / 2 - 



(7) 



where the second identity holds at low frequency. The 
bulk parameter A is given by the usual constraint equa- 
tion, see Chap. 5 of Ref. 9. 

We proceed towards a general analysis of Eqs. (6) at 
the bulk critical point (where the gap V\ cx T vanishes 
as T — > 0). At high temperature, i.e. for T > Tk, we 
have Q = 0. This LM phase is characterized by the Curie 
law xioc(T) = S(S + 1)/T. To determine whether this 
situation persists down to T = 0, we take the Q — > limit 
in Eq. (6), which shows that it does, provided J < J c cx 
e/(S + 1/2) with e = d - 2. The large-A limit therefore 
confirms the existence of a critical interaction for e > 
in the bosonic Kondo model. 

Next we turn our attention to the actual nature of the 
Q 7^ solutions in the regime J > J c and T < Tk- 
Using (7), jl is given from (6) by: 



J 2^ i 



i/0+ 



iv + fi - AQ 2 \v\ d - 2 + 



(8) 



which converges when /i— >0for2<<i<3, and shows 
that there exists a critical value S c of S above which (8) 
cannot be satisfied. Conversely, when S < S ci fi satu- 
rates at T — and cuts the divergence in xioc, signaling 
a complete screening of the local moment (in analogy to 
the fermionic Kondo problem). Therefore, this ES phase 
occurs in the bosonic model for a range of the impu- 
rity spin S < S c , in agreement with our previous strong- 
coupling arguments: the bulk particles are bosonic spin- 
1/2 spinons, which can combine such as to compensate 
higher spins values (up to a maximum value S c ) . Finally, 
when S > S c , fi vanishes and xioc diverges at T = 0, sig- 
naling an undercompensation of the impurity spin. These 
results are summarized in the zero-temperature phase di- 
agrams shown in Fig. 2a. All the previous analysis is well 
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FIG. 3: (color online) Bond parameter Q (lower curve) and 
local spin susceptibility (upper curve) as a function of tem- 
perature for e = 0.5, 5" = 0.3 < S c = 0.5 and J = 2.4. 
Numerical calculations were done with a momentum cutoff 
A = max(fc) = 10 at the bulk critical point g — g c . The 
dashed line denotes the free Curie law \ — S(S + 1)/T. 

borne out by a full numerical solution of the saddle-point 
equations. The existence of an exactly screened phase is 
reflected in the local susceptibility which saturates at low 
temperature, showing a behavior strikingly similar to the 
fermionic Kondo effect, see Fig. 3. 

So far, the large- N discussion applies to space dimen- 
sion d > 2. The case d = 2 is special. First, J c = 0, i.e. 
the local-moment phase is absent (the bulk anomalous 
dimension rj vanishes at N = oo). Second, neglecting 
the spatial dependence of A(x) induced by the impurity 
leads to a loss of the solution of the mean field equations 
at low temperature, an artifact cured through the spatial 
variation of A(x). Although solving for the full A(x) is 
numerically challenging, we have established [by means 
of a systematic small-Q expansion of (6)] that it remains 
short range. Indeed, the effective propagator of the bulk 
bosons, after integrating out the impurity, reads: 

' \ g iu + fi ) 

so that the impurity effect is counter-balanced by choos- 
ing A(x) ~ —(gQ 2 /fi)5(x). We have checked that this 
ansatz leads indeed to a stable numerical solution of the 
mean field equations. The resulting phase diagram is in 
Fig. 2b, and the susceptibility is similar to the one shown 
in Fig. 3. 

Away from bulk criticality. Finally we comment on 
the impurity physics away from the bulk critical point. 
Moving inside the gapped paramagnet, g > g c , shifts the 
critical Kondo coupling J c towards larger values (in par- 
ticular J c becomes non-zero even for e < 0), but screening 
remains at low temperature, an effect which is correctly 
captured by our mean-field calculation. (If one is inter- 
ested in the vicinity of the deconfined QCP described in 
Ref. 14, spinons recombine into triplet magnons above 
the valence bond solid phase, and our large-TV calcula- 
tion does not apply below a typical confinement energy 
scale.) Inside the ordered phase, g < g c at low T, the 
arguments of Ref. 6 apply. 

Conclusions. We have outlined a bosonic version of 



the Kondo effect in magnetic systems with deconfined 
bosonic spinons. Together with the result of Ref. 6, we 
conclude that magnetic impurities are sensitive to the 
nature of excitations in quantum critical magnets, and 
may be employed as a probe of deconfined criticality 
[14]. Applied to CS2CUCI4 we propose to study mag- 
netic properties of dilute impurity moments (e.g. by local 
NMR probes) - true screening would imply the presence 
of bulk spinons. Finally, our results apply to the descrip- 
tion of impurity effects [15] in multicomponent bosonic 
gases trapped in optical lattices [16], albeit with some 
important modifications due to the lack of relativistic in- 
variance and the possible presence of a Bose condensate. 
Interestingly, a quadratic dispersion would imply a bare 
dimension e' = (d — 2)/2 of the Kondo coupling, so that 
a non-trivial critical point could persist in d = 3. 
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